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Abstract 

The two-step contribution to the spin-longitudinal and the spin-transverse cross sections of 
-'^^C,^'^Ca(p,n) reactions at 494 MeV and 346 MeV is calculated. We use a plane-wave approxi- 
mation and evaluate the relative contributions from the one-step and the two-step processes. We 
found that the ratios of the two-step to the one-step processes are larger in the spin-transverse 
cross sections than in the spin-longitudinal ones. Combining these results with the distorted-wave 
impulse approximation (DWIA) results we obtained considerable two-step contributions to the 
spin-longitudinal and the spin-transverse cross sections. The two-step processes are important in 
accounting for the underestimation of the DWIA results for the spin-longitudinal and the spin- 
transverse cross sections. 

PACS nunibers: 24.70. +s,24.10.-i,25.40.Kv 
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I. INTRODUCTION 



In nucleoli induced intermediate energy reactions, a single-step reaction model is the 
simplest way to understand the complex interaction between nucleon and nucleus. However, 
it is insufficient to explain the cross sections at large angles 0, |^. Various groups ^ 
introduced multistep direct processes to compensate for the lack of cross sections in these 
regions. They calculated angular distributions of the two-step and three-step direct 
processes [Q and found that multistep direct processes contribute significantly towards the 
understanding of nucleon induced intermeditate energy reactions. 

The two-step processes have also been applied |^, ^ to the long-standing Rj^/Rt^ prob- 
lem in spin-isospin excitations. With the random phase approximation (RPA), the spin- 
longitudinal response function Rj^ is enhanced and the peak of its energy spectrum is shifted 
downwards, while the spin-transverse response function R^ is quenched and its peak is 
shifted upwards 0]. The theoretical ratio Ri^/R^ should be greater than 1. These quantities 
are extracted experimentally from the measurement of the polarization transfer coefficients 
P, ^, |ll|. With these observables, the ratio is found to be less than 1 and is inconsistent 
with theoretical predictions. 



The spin-longitudinal cross section IDq [0, which corresponds to -Rl |T3|, is roughly 
reproduced by the distorted-wave impulse approximation (DWIA) with the RPA correlation 
in the lower energy transfer region as shown in Fig. |1], but as the energy transfer increases the 



theoretical results become smaller than the experimental results . In the spin-transverse 
cross section I Dp, which corresponds to R^, theoretical estimates reproduce only about half 
of the experimental result in the whole energy transfer region |]14| . 

Nakaoka and Ichimura P] calculated the two-step contribution to IDq and IDp of the 
{p, n) reactions within the framework of the plane- wave approximation to evaluate the rel- 
ative contributions from the one- and the two-step processes. Theoretical results including 
two-step contributions were closer to the experimental results than the DWIA predictions 
for IDq, but they still underestimated IDp. 

In this formalism we assumed that the path lengths of the incident particles in the 
target nucleus for the one-step and the two-step processes are nearly equal, and the effects 
of absorption in the one-step and the two-step processes are similar. Since we wanted to 
evaluate the relative contributions from the one-step and the two-step processes, we removed 
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the effect of absorption from both processes. 

Because of the singularity that appears in the Green's function in the two-step processes, 
numerical integration is difficult. An on-energy shell approximation to the Green's function 
was used as shown in our previous paper However, since this approximation may not be 
reliable it is not used herein. 

In this paper, we utilize the finite e in the denominator of the Green's function in order 
to avoid the singularity and calculate the numerical integration beyond the on-energy shell 
in the transferred momentum space. Then the results are obtained by extrapolation in the 
limit of zero e. In Sec. |T|, we briefly review the formalism of the one-step and the two- 
step cross sections. In Sec. |I|, we present the practical method of numerical calculation. 



In Sec. the results of the one-step and the two-step cross sections in the plane-wave 
calculation are shown. The two-step cross sections with distortion are roughly estimated 
using the DWIA results and they are compared with experimental results. In Sec. 0, we 



analyze the numerical results. Finally, these results are summarized in Sec. VI 



II. FORMALISM 

Consider {p,n) reactions in which the target is excited to the continuum region, for which 
the multiple scattering formalism of Kerman, McManus and Thaler(KMT) may be used. 
In order to see the relative contributions from the one-step and the two-step processes, we 
use a plane-wave approximation. The total Hamiltonian H is written as 

H = Ho + V, (1) 
Ho = + H^, (2) 

where Hp and H^ are the projectile and the target Hamiltonians, respectively, and V is the 
sum of the two-body interactions between the projectile and the target nucleons. 

A. One-step processes 

The one-step T matrix is written in the KMT formalism as 

T«(kf,kO ^ -^(kf|($„|^E^.|$o)|ki), (3) 
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where are the eigenstates of the target nucleus with the excitation energy E^^^^; ki, kf 
are the incident and final momenta, respectively; ti is the transition matrix for the projectile 
and the i th nucleon in the target; and A is the target mass number. 

In the impulse approximation, tj is approximated by the free nucleon-nucleon transition 
matrix (NN t matrix) ^nn- Further in the tp approximation the one-step T matrix is rewritten 
as 

TS{k,,k,) = tNN(q)($n|p(q)|$o), n ^ 0, (4) 

A 

where q = kf — kj is the transferred momentum, and p(q) = ^ e^^^'^^ is the density operator 

i=l 

in the momentum space without the spin dependence. 

To analyze the spin observables we inevitably treat the spin dependence explicitly. We 
introduce the unit vectors 

q . ki X kf „ „ „ 

^ = M' " = ^^' P = ^""- 

In the coordinate system [q, fi, p] the NN t matrices can be decomposed as 

MM 

with fi,p = u, q, n,p, where aju = Ij, o-jfj, = <7j ■ fi {fi u,j = 0,i). The operator ctq and cTj 
denote the spin operators of the projectile and the i th nucleon in the target, respectively. 

The spin-longitudinal cross section IDg and the spin-transverse cross section IDp of the 
{p,n) reactions can be written as [^, |l^ 



ID, = i^^t;(q)i?,,(q,^-*)t,,(q), (7) 
ID, = i^^t;^(q)i?,,(q,u;-*)tpp(q), (8) 
K = (9) 

where Mr is the invariant mass of the residual nucleus, ^/s is the total energy of the system, 
u;™* is the energy transfer between the intrinsic states, and /i; (/if) is the reduced energy 
of the projectile (ejectile). The factor \duj^'^^ / duj\ = y^/M^ is the variable transformation 
coefficient between the energy transfer u in the center of momentum frame and cj™*, and 
i?^/^(q, cu™*) is the response function of the spin density fiuctuation 



i?,,,(q,.;-*) = 5:($o|p^,(q)|«l>„)($„|p,(q)|$o)5(a;-* - i^^), (10) 
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A 

where p^(q) = '^e.'^^'^^aifi. The response functions Rqq and Rpp are nothing but the 

i=l 

spin-longitudinal and the spin-transverse response functions -Rl and i?T, respectively. The 
nuclear correlations are not included in this paper. 



B. Two-step processes 

In the two-step processes the T matrix is written as 



tfo^(kf,kO = -— - E / 7^(kf|($„|— -Eid<^>n')|k„.)xG.,(M(k„.|($n'|^-E^^I*0)|ki). 

(11) 

Since we adopt a plane-wave approximation, the Green's function is diagonal with respect 
to km as 

Gn'{k^) = , , ^ , (12) 

E - y'M2 + - y^(MT + E'^^f + + ze 
where M is the nucleon mass and Mt is the target mass. 

The target gets excited after the first collision, but assuming that the particle-hole pair 
created in the first collision has nothing to do with the second collision, we get IDq and ID^ 
as [|6l 



IDq = kJ duj, J J E E ^i'LKm^'^'(^oK^OM;,^.)^Tr(a,ao^,ao^J,(13) 

respectively, where qi = k^ — ki, uJi is the first-step energy transfer in the center of mass 
frame, and 

- (^)' E E ^k-^^ - - qi)^.^A.(q - qi> q - q'l; - 

X G*(ki + qi; OG(ki + q,; Ot^, (q'Jt,,,, (qi)i?;,;,, (qi, q^; <) (15) 

with /ii, /ii = M, gi, ni,pi and /i2, /i2 = u, q2, n2,p2- Here we define the unit vectors 

q-qi , kf X q2 „ . . , . 

q2 = ] 7, n2 = P2 = q2 X n2. (16) 
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The excitation energy in the Green's function was replaced with the first-step intrinsic energy 
transfer wj^^*, and the response function is extended to a nondiagonal form with respect to 
the transferred momentum q as 

i?,,,(q,q>-*) = 5:($o|p^,(q')|$n)($n|p,(q)|<fo)5(^^°* - E^^'). (17) 

The factor {A — 1)/Ain Eq. ( p!5[ ) represents that the struck nucleon in the first step is never 
struck again in the second step. 

The interference terms with T^g^ and T^g^ are neglected. Exact target states are the sum 
of Ip-lh states, 2p-2h states, and so on. The coefficients of Ip-lh states are expected to 
be random with respect to the exact target states for a fixed Ip-lh state. Then the target 
states can be represented by Ip-lh states with a statistical distribution function. 

It can be considered that the different order T matrices mainly excite different number 
particle- hole states. Then, the final state of T^q"* is different from that of . Hence the 
contribution from the interference terms can be neglected. 

In Fig. ^ the integration region in the transferred momentum space is shown. Outside 
the spheres drawn with solid lines the values of the response functions are negligible. With 
this property of the response functions the integration region is restricted within the overlap 
region of the two spheres. 

In the limit of e = 0, the Green's function is rewritten as 

G (k^;0=P , 

E-^M^ + kl- ^(Mt + ujf')^ + kl 

- in6 {e - + k^ - y^(MT + uT'? + k^) . (18) 

In our previous paper the principal part of the Green's function is neglected (the on-energy 
shell approximation). With the 5 functions the integration region was restricted onto the 
surface of the sphere drawn with a dashed line in Fig. ^. 

Since the validity of this approximation is questionable it is removed. We first calculate 
the sevenfold integration in Eqs. (|13D and (0) for several finite values of e in the Green's 
function to avoid the singularity that appears at e = 0. We then extrapolate the results to 
those at e = 0. 
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III. NUMERICAL CALCULATION 



A. Suppression of energy transfer dependence 

In Ref. [^] the process in which the incident particle loses the energy uji in the first 
collision and propagates in the intermediate state at the energy that is smaller by uJi than 
the incident energy is described, but in this paper the Ui dependence is removed from the 
NN t matrices in the second step and from the Green's functions to make the sevenfold 
integration easier. 

To justify this approximation we show in Fig. ^the Ui dependence of the two-step cross 
sections with the on-energy shell approximation. The removal of the ui dependence from 
the NN t matrices slightly decreases the two-step cross sections for both IDq and IDp. 
Further removal of the ooi dependence from the Green's functions makes them increase a 
little. In spite of these procedures the total variations in the two-step cross sections with 
the on-energy shell approximation are negligible. 

The suppression is also expected to have little effect on the estimation of the two-step cross 
sections without the on-energy shell approximation, because the on-energy shell approxima- 
tion is a restriction of the integration region with respect to the transferred momentum space 
and has nothing to do with the integration of the energy transfer Ui. 

B. Parameters and mesh sizes 

The two-step cross sections are calculated at e = 5.0, 10.0, 15.0, and 20.0 MeV. The 
mesh sizes of the numerical integration are almost the same as those in Ref. 0, although 
the z axis is set parallel to the direction of q. Acu™* = 10.0 MeV, Aq'j^^ = 0.2(l/fm), where 
= {(A - 1)/A}qi, and A cos = 0.125. We set A0i = n/96, 7r/48 for e = 5.0, and 10.0 
MeV, respectively, and A0i = 7r/24 for e = 15.0, and 20.0 MeV. 

IV. RESULTS 

Figure § shows the two-step cross sections of ^^C(j9,n) reaction at 494 MeV as a function 
of e at the energy transfer of 85 MeV. By extrapolating the calculated two-step cross sections 
toward e = 0.0 MeV, they are expected to be 1.90 and 1.93 times as much as the results 
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with the on-energy shell approximation. 

Figure |^ shows the cross sections of ^^C(p,'n,) reaction at 494 MeV as a function of the 
energy transfer in the laboratory frame, ci^iab, where the two-step cross sections without the 
on-energy shell approximation are extrapolated to e = 0.0 MeV. The two-step cross sections 
increase as the energy transfer increases. The results with the on-energy shell approximation 
are shown with dotted-dashed lines. 

The ratios of the two-step cross sections to the one-step ones at c^iab = 125 MeV are 
about 47% and 114% for IDq and IDp, respectively. Those with the on-energy shell ap- 
proximation are revised to include the factor {A — 1)^/A^ in Eq. (^) and they are 24% and 
59%, respectively. The relative contribution is much larger in IDp than in IDq. These data 
are displayed in Table |. 

The cross sections of ^'^C{p,n) reaction at 346 MeV are shown in Fig. ^. The ratios are 
59% and 116% at cjiab = 125 MeV, and those with the on-energy shell approximation are 
34% and 66%. 

The cross sections of ^''Ca(p,n) reactions at 494 MeV and at 346 MeV are shown in Figs. |^ 
and ^ respectively. The ratios at cuiab = 121 MeV are shown in Table |. 

Since the two-step cross sections are calculated in the plane-wave approximation, it is 
difficult to know whether these results account for the underestimation of IDq and IDp. 
One possible way of showing the two-step effects may be to multiply the uncorrelated DWIA 



results with the ratios of the two-step cross sections to the one-step cross sections in the 
plane-wave approximation. Then these renormalized contributions are added to the DWIA 
results with the RPA correlation, since the response functions without the RPA correlation 
were used in the plane-wave calculation of the one-step and the two-step cross sections. 

The two-step contributions for the ^^C(p,n) reaction at 494 MeV thus obtained are shown 
in Fig. ^ The numerical results get closer to the experimental result both in IDq and IDp 
than the previous estimation with the on-energy shell approximation. The two-step cross 
sections have contributions comparable to the discrepancy between the DWIA results and 
the experimental ones. 

The result for the ^°Ca(p,'n,) reaction at 494 MeV is displayed in Fig. |10[ The calculated 
results approach the experimental results both in IDq and IDp and account for a significant 
portion of the previously underestimated spin-transverse cross section. Here one must note 
that the present estimation of the two-step contribution with distortion is crude, but Figs. 



8 



and |T0| indicate to what extent the two-step processes are important to account for the 
discrepancy between the DWIA results and the experimental results. 

V. DISCUSSION 

In this section a qualitative explanation for the origin of the difference between the results 
with and without the on-energy shell approximation is considered by first defining 

^ J-^sJ (^^NN(q - qi)tNN(q - q'i)i?(q - qi, q - q'l; - O 
X G(ki + qOG^ki + q;)tNN(qi)tNN(q'i)^(qi, q'l; (is) 

which is Eq. (|T3|) or (0) but the integration with Ui is not included. Here we suppressed the 
spin dependence for simplicity. Because the on-energy shell approximation is a restriction 
of the integration region in the transferred momentum space, the accuracy of the on-energy 
shell approximation is considered only with regards to the integrations of the transferred 
momenta. 

The quantity X^"^^ has to be the square of the absolute value of a complex number, 
because the cross section is obtained from the square of the T matrix. Actually the response 
function is the sum of target excited states as written in Eq. ([T7| ) and each term is divided 
into a function of qi and a function of q'j^, and Eq. ( |T9|) is rewritten as the product of the 
integration of qi and its complex conjugate, i.e., the integration of q'^^. 

Dividing the Green's functions into the real part and the imaginary part, one gets 

^^'HO = / / (^^NN(q - qi)tNN(q - q'i)^(q - qi, q - q'l; ^'"^^ - O 

X {ReG(ki + qi)ReG'(ki + q^) + iImG(ki + qi)ReG(ki + q'^) 
- iReG(ki + qi)ImG(ki + q'^) + ImG(ki + qi)ImG(ki + q^)} 

X tNN(qi)iNN(q'i)^(qi,q'i;<') (20) 

It is the fourth term that remains when the on-energy shell approximation is applied. The 
same argument as in the previous paragraph can be made with respect to the first term and 
the fourth term, and both of these two terms become positive. The second and the third 
terms are the complex conjugate of each other, and the sum gives a real number. However, 
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TABLE I: The ratios of the two-step cross sections to the one-step ones at 125 MeV for ^^C and 
at 121 MeV for ^''Ca. Values in parentheses are those with the on-energy shell approximation. 

Target iCiab(MeV) IDg{%) IDp{%) 
494 47(24) 114(59) 

346 59(34) 116(66) 

40Ca 494 96(46) 192(88) 
346 92(50) 172(89) 

the contributions to IDg and IDp from the second and the third terms are found to be small. 
In Fig. I, they are only -1.7 x 10"^ and -2.7 x 10"^ out of 1.29 x 10"^ and 1.08 x 10^^ 



mb at e = 20.0 MeV, respectively. Therefore the cross sections obtained in Sec. |^ are larger 
than those with the on-energy shell approximation. 

The contributions from the fourth term are displayed in Fig. ^ with dotted lines. Their 
extrapolation toward e = 0.0 MeV seem to hit the results with the on-energy shell approx- 



imation. This shows that the results obtained in Sec. |^ are consistent with the previous 
estimations with the on-energy shell approximation. 



VI. SUMMARY 



Two-step contribution to the intermediate energy {p,n) reactions were calculated using the 
plane-wave approximation. We carried out the integration beyond the on-energy shell. The 
ratios of two-step to one-step contributions were found to be larger for the spin-transverse 
cross sections than for the spin-longitudinal cross sections. The two-step cross sections 
without the on-energy shell approximation are about twice as large as our previous results 
with the on-energy shell approximation. 

The cross sections including two-step processes were compared with experimental results. 
We conclude that the two-step cross sections have contributions comparable to the discrep- 
ancy between the DWIA and the experimental results in the highly excited region for IDg 
and they account for a portion of the previously underestimated IDp. 
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FIG. 1: The DWIA results for the ^'^C{p,n) reaction at 494 MeV |T|]. The upper panel shows the 
spin-longitudinal cross section, and the lower panel shows the spin-transverse one. The horizontal 
axis is the energy transfer in the laboratory frame. 
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FIG. 2: The integration region in the transferred momentum space. Two spheres drawn with soUd 
Hncs have a radius 4.8(l/fm) for ^^C, beyond which the response functions are negligible. The 
dashed sphere is the on-energy shell. 
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FIG. 3: The dependence of the first-step energy transfer u>i for ^^C{p,n) reaction at 494 MeV. The 
dotted fines indicate the two-step cross sections with the fufi oJi dependence. The dotted-dashed 

fines are those without the uji dependence in NN t matrices, and the dotted-dotted-dashed fines 
represent the two-step cross sections without the ui dependence both in the NN t matrices and 
the Green's function, respectively. 
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FIG. 4: The two-step cross sections of ^^C(p,n) reaction at 494 MeV as a function of e. The energy 
transfer is 85 MeV. The dashed lines indicate the results extrapolated smoothly from the calculated 
results. The dotted lines do the contributions from the fourth term in Eq. (pOl). The results with 
the on-energy shell approximation are plotted with cross symbols at e = 0.0 MeV. 
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FIG. 5: The cross sections of ^^C(p,n) reaction at 494 MeV as a function of wiab- The two-step 
cross sections are indicated with dotted hnes, and they are extrapolated to e = 0.0 MeV. The 
dotted-dashed lines are those with the on-energy shell approximation. The dashed lines are the 
one-step cross sections. 
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FIG. 6: Same as Fig. |, but at 346 MeV. 
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FIG. 7: Same as Fig. |, but for ^^Ca. 
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FIG. 8: Same as Fig. |, but at 346 MeV. 
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FIG. 9: The DWIA with the two-step results for the ^^C{p,n) reaction at 494 MeV. The upper 
panel shows IDq, and the lower panel does I Dp. The dashed lines indicate the DWIA results. The 



solid lines include the two-step effects. Experimental data are taken from Ref. |10| 
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FIG. 10: Same as Fig. |, but for ^OQa. 
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